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Abstract 

We calculate heavy quark symmetry breaking in the slopes and curvatures of the B D^*^ii> 
spectra at zero recoil, including the order a^Po corrections. We point out that the theoretical 
uncertainties in the differences between B ^ D and B D* slopes and curvatures are smaller 
than in the deviations of the slopes and curvatures themselves from their infinite mass limits. We 
find that the central values of the current experimental results for the difference of the slopes differ 
from our calculations when QCD sum rules are used to estimate subleading Isgur-Wise functions. 
A better understanding of the shapes of the B D^*'>li> spectra may also help to reduce the error 
of \Vcb\ extracted from the zero recoil limit of i? — > D*iv. We argue that heavy quark symmetry 
requires that the same fitting procedure be used in the experimental determinations of the shape 
parameters and \Vcb\ from the B DlU and B D*(.v spectra. 
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I. INTRODUCTION 



The determination o 
heavy quark symmetry 



" \Vcb\ from exclusive B D^*^^ decays is based on the fact that 
relates the form factors which occur in these decays to the Isgur- 
Wise function, whose value is known at zero recoil in the infinite mass limit. The symmetry 
breaking corrections can be organized in a simultaneous expansion in ag and Aqcd/'^q 
{Q = c,b). The B — > D'^*Hi' decay rates are given by 



487r3 



■r3(l-r,)Vw2-l(w+l)^ 



4w 1 — 2wr* + rl 

1 + 



1 + w (1 - r*)2 

:2 ^5 

dw 487r3 

where w = v ■ v' = {m\ + — q'^)/{2mBm^(,)) and r(^.) = m^(,)/mB- Both J^(ty) 

and J-'^{w) are equal to the Isgur-Wise function in the mg — ^ oo limit, and in particular 
= JF^(I) = 1, allowing for a model independent determination of \Vcb\- 

The main theoretical uncertainties in such a determination of \Vcb\ come from the value 
of and from the shape of T(^^){w) used to fit the data. Such a fit will continue to 

be important, since the number of -B ^ D*^z/ events needed to measure |Vc6|J^*(l) with 
a statistical error of order (Aqcd/'^q)^ scales parametrically as (mg/AQCD)^- This can 
be seen from Eq. (0), making no assumption on the shape of J-'^{w), by considering a bin 
at zero recoil with width of order of the desired accuracy, that is, of order (Aqcd/'^q)^- 
Similarly, when unquenched lattice QCD calculations of will be available with error 

a, the number of events needed to measure \Vcb\ with a comparable statistical error, 

without assumptions about the shapes of will scale as a"^/^ in i? ^ D*^z/ and 

as a"^/^ in i? ^ Diu. Reliable unquenched lattice QCD results for JF(^-)(1) are likely to 
be available before comparable results for the functional form of Constraining the 

shapes of ^(^^^{w) will remain important in the foreseeable future. 

The zero recoil limit of J-'(^^){w), including symmetry breaking corrections, can be written 
schematically as 

jc-^(l) = 1 + CA{as) + — + ^ + . . . , 

rriQ 

^(1) = l + cv(a.) + ^ + y + .... (2) 



Theperturbative corrections, ca = —0.04 and cy = 0.02, have been computed to order 
a1 j^, and the unknown higher order corrections should be below the 1% level. The order 
^qcd/'^q correction to JF^(l) vanishes due to Luke's theorem The terms indicated by 
(. . .) in Eqs. (0) are only known using phenomenological models or quenched lattice QCD at 
present. This is why the determination of \Vch\ from B D*£ij is theoretically more reliable 
for now than that from B Diu, although both QCD sum rules and quenched lattice 
QCD 0] suggest that the order Aqcd/^^q correction to is small. Due to the extra 

— 1 helicity suppression near zero recoil, B — > Diu is also more difficult experimentally. 

Some of the order AqcD/fng corrections which enter jF(^)(w) also influence ratios of form 
factors measurable in B ^ D*iu decay. The exclusive semileptonic B — > D*ii' decay rate 
is parameterized by four form factors. 



, hA, {w + 1) e*^ - {hA.v^^ + hA.v'^') (e* ■ v) . (3) 

y/mo* rriB 

One linear combination is not measurable when the lepton masses are neglected. The form 
factor hAi dominates the rate near zero recoil. It is conventional to define two measurable 
ratios of the form factors 

hyjw) hA;{w) + {mD,/mB)hA^{w) 
= - — -— , K2{w) = — . (,4j 

hAAw) riAAw) 

In the infinite mass limit Ri{w) = R2{w) = 1, and deviations from this limit measure certain 
combinations of the subleading Isgur-Wise functions as it will be discussed it later. 

Since there are several form factors in i? ^ D*iu, it is customary to fit the shape 
parameters of /i^i, together with the form factor ratios i?i_2. For the shapes of J-'{w), J-'^{w), 
and analyticity imposes stringent constraints between the slopes and curvatures at 

zero recoil LZll. It is convenient to write 



J^{w) = ^(1) [l -p%{w-l) + ct - 1)' + ■ ■ 
T,{w) = [l - p^^ {w-l) + c^, {w - 1)2 + 

hA, (w) = hA,{l) [l - {W-1)+ CA, {w-iy + ---\. (5) 

A lower index X will denote any of JF, JF^ or Ai. In the rriQ oo limit, heavy quark 
symmetry predicts px = Po and cx = Cq, the slope and curvature of the Isgur-Wise function, 
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FIG. 1: The measured values of |V^fe|.F*(l) and p\_^ from unitarity constrained fits to ^ D*£u 
decay (left), and \ Vcb\ ^(1) and pjr from linear fits to the B Dii? decay (right). From Ref. Q]. 

respectively. These equalities are violated by corrections at order and Aqcd/''^q- In a 
linear fit to the data the expansions in are terminated at linear order in w — 1, while in a 
quadratic (or "free curvature") fit the expansions are terminated at order {w — 1)^. Higher 
order terms are small but non-negligible, since the fitted range extends to w = 1.59 (1.50) 
in 5 — s> D^*^ decay. Unitarity constraints give a strong correlation between pjr^ ^ and the 
coefficients of the higher order terms in w — 1 0? and also constrain the effect of the 
higher than quadratic terms in w — 1 . For these reasons, an unconstrained quadratic fit 
to the data is not equivalent to a unitarity constrained fit. 

In the next section we compute the violations to the symmetry relations between the 
slope and curvature parameters and find them to be small. Moreover, we observe that 
in the present data the slope parameters pjr^ and differ significantly if the spectra are 
fit linearly or quadratically; see Table I and Fig. 1. However, when present data is fit to 
the functional form constrained by unitarity, the result is consistent with the symmetry 
relation = p'jr. Therefore, if in the future the data supports the necessity of a unitarity 
constrained form for jF^(w), then it will be equally necessary to use a unitarity constrained 
fit for T{w). 

Furthermore, to the extent that we can reliably estimate the deviations from pjp^ = p^r. 
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Fitted slope parameter 


CLEO 


BELLE 


B - 


D*l9, unitarity constrained fit to 


1.67 ±0.11 ±0.22 [10] 


1.35 ±0.17 ±0.19 [nj 


B - 


D*l9, linear fit to pjp^ 


0.98 ± 0.09 ± 0.07 [12] 


0.89 ±0.09 ±0.05 [13] 


B - 


-> D(.9, unitarity constrained fit to pjp 


1.30 ±0.27 ±0.14 [14] 


1.16 ±0.25 ±0.15 [15] 


B - 


Dip, linear fit to p'jr 


0.76 ±0.16 ±0.08 [14] 


0.69 ±0.14 ±0.09 []J 



TABLE L The most recent available CLEO and BELLE measurements of the slope parameters, 
Px, in i? — > D^*^£i> decay. The difference between fitting to the unitarity constrained shapes of 
one or the other of Refs. |7| or is very small. 



the accuracy in the extraction of \Vcb\ can be improved by a simultaneous fit to i? ^ D* and 
B ^ D. While the value of \Vch\ -^(1) has larger uncertainty than \Vcb\ since B — * Dip 

is helicity suppressed near zero recoil, the error on ^ is comparable as extracted from 
B D* and B ^ D. Since there is a very strong correlation between ^ and \Vcb\ J^{*)(1), 
despite the fact that J^(l) is less well-known than J?^*(l), a better determined slope may 
have significant implications for the value of \Vcb\ as it is extracted from the zero recoil limit 
of -B — s> D*ih'. When precise unquenched lattice calculations of J^(^){w) become available, 
consistency between the predictions and data for px niay be an important cross-check of 
the \Vcb\ determination. 



II. ANALYTIC RESULTS 

We next give the theoretical predictions for the heavy quark symmetry violation in p^, 
cx, and the form factor ratios R12 at order Aqqd/^q, C(s and a^Po. The first two corrections 
are known in the literature, but the order a^/3o terms, which probably dominate the order 

corrections (since (3o = 11 — 2?2//3 is large) are new. These are required to predict heavy 
quark symmetry breaking in the quantities under consideration at the ~ 5% level, and may 
become important in testing lattice results. We can write the slope parameters p|^ as 

Px = pI + ^'x'' + ^/x'"''^ (6) 

where X = JF, JF^, Ai denotes the functions under consideration. The perturbative correc- 
tions can be calculated model independently and are contained in ^x'"*, while 6x^"^^ contains 
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the order Aqcd/"^q corrections to p\. The perturbative corrections can be computed in 
eluding the order a^/^o corrections from expanding Eqs. (39) in Ref. 



3. 



We obtain 



^JF, — 



as 

TT 



In z 



9(1-^)2 9(1 
2 rC_17^ + 3^2 (i + ;2)(i2-23z+12z2) 



27(1-^)2 



54(1 



^3 



In z 



_ 5K) 



4(1-^ + ^2 
9(1 -z)2 



+ 



22(1 + 2) 



9(1 



In 2; 



4 + 52 + 4^2 (1 + 2) (12 - 37z + 12^2 



In z 



(7) 



54(1 -z)2 108(1-^)3 

where z = mc/mi), and denotes the strong couphng renormalized in the MS scheme at the 
scale /i = y/mf,mc. The terms of order agree with Ref. |l7|. Using z = 0.29, cXg = 0.26, 
and Po = 25/3, we obtain 



6 



0.079 + 0.046 



6 



K) 



(5 



K) 



0.034 + 0.018, 



(8) 



where the first and second terms come from the order and a'^Po corrections, respectively. 
The apparent bad convergence of the perturbation series in Eq. (jH)) may be due to the 
fact that they contain so-called renormalon ambiguities, and are only well-defined physical 
quantities when the nonperturbative corrections discussed next [see Eq. (P)] are included. 
This assertion is supported by the fact that we will find significantly better behavior when 
these series are expressed in terms of a short distance mass in Sec. III. 

The corrections in 6x^^^ depend on the four subleading Isgur-Wise functions that pa- 



rameterize first order deviations from the infinite mass limit. Using the notation of 
finds from Ref. 



one 



r{l/m) 



5(1 + 2) 



16 



6 

1 + z 



X2(l) - 16x'3(l) + 



2z 



3(1 - z 



5z^ ,..2(1 
V{1) + 



l + z 



^4 , 1 + 2: + 2^2 

+ 3^^^^^+ 3(1-.) ^^^^ 



(9) 



where prime denotes d/dw. rj{w) = iz{w)/^{w) parameterizes the (9(AQCD/"^c,b) correc- 
tions to the b c current, and X2,3{'w) describe the matrix elements of the time ordered 
product of the chromomagnetic operator, {Qs/'^) h^cr^uG^'^h^, with the leading order current. 
An important point is that the poorly known function Xi{'^) which parameterizes matrix 
elements involving the time ordered product of the kinetic energy operator, hv{iD)'^hy, with 
the leading order current drops out from the differences in Eq. (jU)). In general, Xi does not 



effect any quantity determined by heavy quark symmetry at leading order, such as the slope 
and curvature differences in Eqs. and (fT^ . and i?i,2 in Eqs. flT^ and (fT^ . The reason 
is that Xi enters all form factors in the combination C,{w) + {A/rric + A/nif,) Xi{w), where 
C,{w) is the Isgur-Wise function. 

It is straightforward to compute the heavy quark symmetry breaking corrections between 
ct, C:f^, and ca^- Similar to Eq. we define 

A 



cx = co + A^x'^ + 



2mr 



A 



(l/m) 
X 



(10) 



For the differences of A^^"^ we obtain 



X 



A 



2as 
1357r 

+ 



47 - 1482 + 282^2 - 148^3 + A7z^ 5z(l + z)(l + 6z + z'^) ^ ' 

\ ^ In z 



4057r' 



509 - 1241z + 1084^2 _ 1241^3 ^ 509^4 
5(1 - z)* 

(1 + z)(282 - 1277^ + 2142^2 _ 1277^3 + 282z^) 



A?;^ - A^^ 



2a. 



277r 
+ 



4(1 - zf 

{as) Aas)' 

7-23Z + 12z^ - 23^3 ^7^4 z{l + z){l - I2z + z^) 



\nz 



32477" 



(1 - zY 
44 - 83^ - 134^2 



+ 



^3^3 + 44^^ 



l-z) 



In 2; 



{i-zY 

[I + z)(42 - 187z + 396^2 - I87z=^ + 42z^) 



4(1 -z) 



Inz 



- 



(11) 



The order as piece in the first of these equations agrees with Ref. [18|. To this order, pi can 
be taken as any of p^, p^^, or p\^. With the previously used values of z and a^, we obtain 



= (0.074 + 0.043) + pI (0.079 + 0.046) 
A^";^-A5?;^ = (0.058 + 0.031) + p^ (0.034 + 0.018) 



(12) 



where the first and second terms in each parenthesis come from the order and a^/^o 
corrections, respectively. The convergence of these series is again quite poor. The order 
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Aqcd/'"^c,6 heavy quark symmetry breaking in the curvature differences is given by 



A 



(l/m) 



A 



(l/m) 



(l + ^)(25 - 42^ + 25^2 



36(1 - zf 



4X2(1) 



1-62 + ^2 



16 



X2(l) + 8x'3'(l) 



+ ^(1) 



5 - 28z + 18^2 - 52^3 + 25z^ 



18(1 -z) 



V(l) 



1 - 2^ + 
3(1-^) 



(1 + ^) 



+ Po 



(l/m) 



16 



X2(l) + 16X3(1) 



, (l/r") 



2(l + 2)(2-3z + 22:^) , ^ n^i^.^i±f! 4 
9(1^^)^ + ^^^^'^ 9(1-.)^ - 3 



+ ^(1) 



5 - 19z - 9^2 - 25z^ 



18(1 - z) 



16z^ ,,^,1 + ^ + 2^2 

^1^71 ^ 

3(1-2;) 



+ Po 



(l/m) 



4 



X2(l) 



(13) 



Note that the coefficients of pg are independent of Xi- While not all parameters entering 
these formulae are known, it may be possible to get information on X2(l) ^"(1) from 
precise measurements of J^^iX) and -Ri(l), respectively, or from lattice QCD calculations. 

Independent of model calculations, experimental data on the B —>■ D*£u form factor 
ratios defined in Eq. ^ will constrain some of the subleading Isgur-Wise functions entering 
Eqs. IQ and (fT^ . Measurements of Ri^2 can be used to constrain the quantities //(I), ''^'(l), 
and X2(l) according to 



^i(l) 
R[{1) 



Atts at ^ 
l + ^ + ^/?o 



2 _ l + z 
9 ~ 3(1 - z) 



hiz 



A 

2mr 



l + z-2zri{l) 



+ ... 



4a, a, 


97r TT 



A 

2mr 



2 

s 

2 

1+Z 



2(1 + ^2) ^ (1 + _4^ + ^2) 

9(1 -z)2 9(1-^)3 
- zri{l) + 2zri'{l)\ + ... . 



\nz 



(14) 



Yet again, one encounters badly behaved perturbation series, 0.110 + 0.055 for -Ri(l) and 
-0.037 - 0.025 for R[{1). For R2 we obtain 



^2(1) 



1 - 



2a, 
37r 
A 



2z zil + z) , ■ 
+ -7^ 777 Inz 



1-z (l-z) 



13a^ 
187r2 



Po 



z z{l + z) ^ ■ 
+ ttt: 777 Inz 



l-z 2{l-zy 



2mr 



(l + 3^)r/(l) + 4(l-z) X2(l) 



+ 



^2(1) 



2a., 



+ 
+ 



37r 

187r2 
A 
2mr 



z{l + lQz + z^) 2z\l + z) , ■ 
H 7z 71- Inz 



Po 



3(1-^)3 (l-z)^ 
z(l- 38^ + ^2) z{l + z){l-llz + z^) 



(15) 



+ 



2(1-^)3 ' (l-z)^ 
'r^(l) 



In 2; 



(1 + 3^) ^ - V(l) - 4(1 - z) (x',(l) + X2(l) pI) 



+ 
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Note that the 0{as) corrections to -R2(l) and -R2(l) very small, around +0.0056 and 
—0.0011, respectively, and the C(a^/9o) corrections are even smaller, +0.0021 and —0.0006, 
respectively. This is not unexpected, since all the nonperturbative corrections to R2 involve 
the subleading Isgur-Wise functions but not A by itself, so the perturbation series does not 
involve a leading renormalon. 



III. DISCUSSION AND CONCLUSIONS 



To evaluate the results of the previous section, they must be expressed in terms of short 
distance quark masses. We use the upsilon expansion and express rric through mj, — 
rric = TriB — mo + Ai/2m5 — Ai/2m/), where = (m^ + 3mB.)/4 = 5.313 GeV and 
rnD = {mr, + 3m£).)/4 = 1.973 GeV. It is convenient to re-express as in terms of as{mh). 
For the IS" h quark mass we use = 4.75 ± 0.07 GeV, which follows from Ref. [2^ using 
the CLEO measurement of {E.y) = 2.346 + 0.034 GeV in B ^ Xg'-f corresponding to a photon 
energy cut > 2 GeV [2I]. The uncertainty in , and the errors related to it which we 
quote below, will almost certainly be significantly reduced in the near future. In the upsilon 
expansion 

K = mB-ml^ - 0.051 e - 0.091 c^lm + • • • , (16) 

where we used asirrih) = 0.22, e = 1 is the parameter of the upsilon expansion, and e^Lj^ 
denotes the part of the second order correction proportional to /?o- Since we include the 
a^Po terms in all results, the residual scale dependence is very small. 

The right-hand sides of Eqs. and ()13p can only be estimated at present using model 
predictions. The subleading Isgur-Wise functions have been computed including order as 
corrections in QCD sum rules, yielding 

X2(l)^-0.04, X3(l)-0.02, r/(l)^0.6, V(1)~0. (17) 

Using the upsilon expansion to eliminate the quark masses, and these values for the numerical 
estimates, we obtain 

p^-p^^ = 0.203 + 0.053 e- 0.013 e|LM + 0.075 77(1) + 0.14^(1) 

+ 1.0 X2(l) - 3.0 X3(l) - 0.018 Ai/GeV^ ^ 0.19 , (18) 
- p^^ = 0.081 + 0.024 e - 0.006 e^LM + 0.131r/(l) + 0.25^2(1) - 0.007 Ai/GeV^ ~ 0.17 , 



We also used Ai = —0.25 GeV^, although the results are hardly sensitive to the value of this 
parameter. The behavior of these perturbation series are clearly much better than those in 
Eq. (jSI) in terms of the quark pole masses. The uncertainty due to a ±70 MeV change in 
ml'^ is =f0.020 and =1=0.022 in these estimates of — pjr^ and p\_^ — p^^, respectively. 

Although these estimates are model dependent, they are less so than one might at first 
think, since the first terms on the right-hand sides of Eqs. 0, which are model independent, 
contribute a large part of the result. These results mostly depend on the value of 77(1) and 
on the smallness of the functions X2,3{'w), which parameterize order Aqcd/^^q corrections 
due to the chromomagnetic operator. If X2(l) and x'^i^) were order unity then these results 
could be dramatically different. However, X2,3{w) are expected to be small in most models 
(recall that X3(l) = due to Luke's theorem jS]), and as we will see below, //(I) can be 
constrained from i?i^2(l)- 

For the curvature differences we obtain 

- = 0.202 + 0.050 e - 0.012 e^LM - 0-087 r/(l) - 0.08 r/'(l) - 0.07 r/"(l) 
-0.12x2(1) -1.0x2(1) + 1-5 X3(l)- 0.011 Ai/GeV^ 
+ Po [P^ - Pk - 1-0 X2(l) + 3.0 X3(l)] ^ 0.17 + 0.29 pi , 
ca,-c^, = 0. 141 + 0.042 e- 0.007 e^LM- 0.059 77(1) -0.13^(1) 
- 0.12 X2(l) - 0.25 X2(l) - 0.005 Ai/GeV^ 

+ Po [Pm - pk - 0.25 X2(l)] ^ 0.14 + 0.18 p2 , (19) 

where for the numerical estimates we also used X2(l) — 0.03 ^, and x'si^) = v"{^) = 0. 
The uncertainty due to a ±70 MeV change in ml^ is =f0.021 and =f0.018 in the 0.17 and 
0.14 terms in these estimates of cjc- — cjr^ and cai — cjc-^, respectively. While there is a 
sizable uncertainty again due to the subleading Isgur-Wise functions, for their particular 
values predicted by QCD sum rules, the final result is dominated by terms which are model 
independent. 

For the form factor ratios Ri and R2 we obtain 

Ri{l) = 1.243 ± 0.079 e- 0.016 e^LM- 0.112 r7(l)- 0.021 Ai/GeV^^ 1.25, 

R[(l) = -0.122 - 0.021e ± O.OlOe^L^ + 0.0567^(1) - 0.112r/'(l) ± O.OllAi/GeV^ ~ -0.10 , 

i?2(l) = 1 ± 0.006 e ± 0.001 e^LM - 0.355 77(1) - 0.53 X2(l) ^ 0.81 , (20) 

R'^(l) = -0.001 e ± 0.178 ri{l) - 0.355 ^(1) - 0.53 Xall) - 0.53 X2(l) Po - 0.09 ± 0.02 pi . 
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The uncertainty due to a ±70 MeV change in ml^ is =f0.03, ±0.01, ±0.03, and =f0.01 in these 
estimates of -Ri(l), -^'1(1)5 -^2(1), and ^^2(1)5 respectively. Unfortunately the sensitivity to 
the values of the subleading Isgur-Wise functions is not too large, since they enter with small 
coefficients. Still, useful constraints can be obtained, e.g., -R2(l) measures r]{l) assuming 
that X2(l) is small (or a linear combination of them otherwise), ri'{l) can be measured using 

Ri{l) + 2R[{1) = 1 ± 0.037 e± 0.004 e^LM- 0.223 r/'(l), (21) 

and then R2{^) can be used to constrain a linear combination of ^2(1) and X2i^)- The 
coefficient of ^^'(1) changes by ±0.032 under a ±70 MeV variation of ml^ , while the other 
terms are essentially unaffected. Since -R'/2(l) are expected to be even smaller than R[ 2(1); 
it seems unlikely that they could give useful independent information. CLEO measured 



Ri = 1.18±0.30±0.12 and R2 = 0.71 ±0.22±0.07 [23, assuming that they are independent 
of w and that hAi{w) has a linear w-dependence. These results agree well with Eq. (EUj) . 

Comparing our results for heavy quark symmetry breaking in the slope parameters in 
Eq. ()18|) to the data in Table I, there are several points to be made: 

(i) The result of the linear fits at both CLEO and BELLE indicate — p^r ~ 0.2 (see 
Table I). This is opposite to what is expected based on the QCD sum rule predictions for the 
subleading Isgur-Wise function in Eq. (fTSjl . The simplest way to accommodate the central 
value of the data is if x'si^) ~ 0.15, which is several times larger than the QCD sum rule 
prediction. It should be straightforward to decide using lattice QCD whether this large value 
of x'si^) occurs. 

(ii) The result of the unitarity constrained quadratic fits at both CLEO and BELLE 
indicate p\_^ — p^ ~ 0.3 (see Table I), whereas based on the QCD sum rule predictions for 
the subleading Isgur-Wise functions one would expect this difference to be close to zero, 

Pm -Pr = -0-14 ± 0.06 r]{l) - 0.14 ^(1) - 0.75 ^2(1) ± 3.0 ^3(1) ± . . . ^ -0.02 . (22) 

The value X3(l) ~ 0-15 suggested above also accommodates this data well. 

(iii) The value of X3(l) is hard to constrain from data, since its contribution to R'i2{^) 
is suppressed by as/vr. A large value of x'zi^) could explain a large heavy quark symmetry 
breaking in the slope parameters, so determining x'zi^) from the lattice is very important. 
Of course, it would be desirable to compute all subleading Isgur-Wise functions from the 
lattice. 
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(iv) Since heavy quark symmetry breaking would be unacceptably large in the comparison 
of obtained from the unitarity constrained quadratic fit with obtained from the 
linear fit, one must use the same fitting procedure to extract the slopes and compare them. 
Therefore, one must be careful not to draw wrong conclusions when comparing the two sides 
of Fig. 1.^ We expect the data for B D*iu will eventually favor the unitarity constrained 
fit, at which point the data for B D£u will have to be fit with the same method lest 
one gives up the symmetry relations that follow from the heavy quark expansion or invokes 
surprisingly large values for subleading Isgur-Wise functions (making convergence of the 
expansion questionable). 

(v) Measurements of heavy quark symmetry breaking in the slope and curvature param- 
eters, p'j^ and Cx, together with measurements of the Ri^2 form factor ratios will strongly 
constrain the order Aqcd/'^q corrections. Eqs. (fTHj) . (fT^ . and ^lU^ can be used to test 
future lattice calculations, or model predictions, such as those from QCD sum rules used in 
this paper for some numerical estimates. 

(vi) Most importantly, a better knowledge of the slope parameters will help to reduce the 
error of |V"cfe|, since there is a very strong correlation, as can be seen from Fig. 1. 

In conclusion, we calculated heavy quark symmetry breaking in the slopes and curvatures 
of the B D^*^^ spectra at zero recoil, including the order a'^Po corrections. A combined 
fit to the shapes of the B D*ijj and B Diu spectra together with the form factor 
ratios -Ri and R2 may lead to a better knowledge of the values of the subleading Isgur- 
Wise functions. This in turn may reduce the error of |Vc6|, since the correlation between 
the p^^^j and \Vcb\ J-'[*){1) is very large. Once p^ — p^^ is better understood, it will also be 
interesting to see how well the cjr — cjr^ constraint is satisfied by the data fitted using the 
unitarity constraints. When jF(^)(w) is computed in unquenched lattice QCD, agreement of 
the predicted values of p^r — p^r^ and cjc- — cjc-^ with data would give additional confidence 
(beyond checking that \Vcb\ extracted from B ^ D and B D* agree) that the errors 
are well understood. It will be especially reassuring if such a value of \Vcb\ agrees with the 
inclusive determination at the few percent level. 

^ As the present authors first did. 
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